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Abstract 

Background: Ebola viruses cause viral hemorrhagic fever in humans and non-human primates and are endemic in 
Africa. Reston ebolavirus (REBOV) has caused several epizootics in cynomolgus monkeys {Macaca fascicularis) but is 
not associated with any human disease. In late 2008, REBOV infections were identified in swine for the first time in 
the Philippines. 

Methods: A total of 215 swine sera collected at two REBOV-affected farms in 2008, in Pangasinan and Bulacan, 
were tested for the presence of REBOV-specific antibodies using multiple serodiagnosis systems. A total of 98 swine 
sera collected in a non-epizootic region, Tarlac, were also tested to clarify the prevalence of REBOV infection in the 
general swine population in the Philippines. 

Results: Some 70 % of swine sera at the affected farms were positive for REBOV antibodies in the multiple 
serodiagnosis systems. On the other hand, none of the swine sera collected in Tarlac showed positive reactions in 
any of the diagnosis systems. 

Conclusions: The high prevalence of REBOV infection in swine in the affected farms in 2008 suggests that swine is 
susceptible for REBOV infection. The multiple serological assays used in the study are thought to be useful for 
future surveillance of REOBV infection in swine in the Philippines. 
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Background 

The Ebola virus (EBOV) is an enveloped, negative- 
strand RNA virus belonging to the family filoviridae in 
the order of mononegavirale [1]. Four of the five ebola- 
virus species, Zaire (ZEBOV), Sudan, Tai Forest, and 
the recently discovered Bundibugyo ebolavirus, are 
endemic in continental Africa and cause a severe form of 
viral hemorrhagic fever with high mortality in humans 
and non-human primates [2-6]. Reston ebolavirus 
(REBOV) is sporadic in the Philippines and has caused 
several epizootics in cynomolgus macaques [7]. REBOV 
was first isolated in 1989 from cynomolgus macaques 
imported from the Philippines for medical research in 
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the United States [7-10]. About 1,000 monkeys died or 
were euthanized in a quarantine facility in Reston, Vir- 
ginia. Subsequently, 21 animal handlers at the Philip- 
pine exporter and four employees of the quarantine 
facility were found to have antibodies to the virus, indi- 
cating that they had been infected [11,12]. Epizootics in 
monkeys in the Philippines were then reported in 1992 
and 1996, and all the epizootics have been traced back 
to a single monkey facility, in Calamba, Laguna in the 
Philippines [11,13-16]. Since the closure of the facility 
in 1997, no REBOV epizootics in cynomolgus monkeys 
have been reported. 

In October 2008, REBOV infection was confirmed for 
the first time in swine associated with multiple epizoo- 
tics of respiratory and abortion-related diseases in the 
Philippines [17]. In several pools of swine samples col- 
lected from geographically distant swine farms, co- 
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infection with REBOV and porcine reproductive and re- 
spiratory syndrome virus (PRRSV) was confirmed [17]. 
Serological studies of limited scale on 13 swine sera in 
the affected farms failed to detect REBOV antibodies in 
ELISA, although PRRSV antibodies were detected [17]. 
It is still unclear how REBOV was spread among swine 
during the epizootic. Moreover, it is not clear if REBOV 
infection in the swine population is either sporadic and 
incidental or common in the Philippines. To try to an- 
swer these questions, we prepared multiple serodiagnosis 
systems for detecting REBOV infection in swine and 
analyzed swine sera obtained from the affected farms 
and from farms not associated with any epizootics in the 
Philippines. The results showed a high prevalence of 
REBOV infection in swine in the affected farms at the 
epizootics in 2008; however, REBOV antibodies were not 
detected in the swine population not associated with the 
epizootics, indicating that REBOV infection in swine in 
the Philippines is not common, at least in some parts 
of Tarlac. 



Results 

Detection of REBOV-NP and -GP antibodies in swine using 
IFA 

Swine sera were analyzed for the presence of REBOV- 
NP and -GP antibodies in IFA specific to REBOV-NP 
and -GP, respectively. In the IFA, none of the 49 swine 
sera collected in Japan showed a positive reaction (data 
not shown), and so they were considered to be 
REBOV-NP and -GP antibody negative. In the IFA spe- 
cific to REBOV-NP, antibody positive swine sera showed 
characteristic granular staining patterns in the cytoplasm 
(Figure 1A), which were indistinguishable from those of 
REBOV-infected cynomolgus monkey sera [18] and 
REBOV-NP immunized rabbit sera (data not shown). 
Antibody-negative swine sera showed no reaction (Fig- 
ure IB). In the IFA specific to REBOV-GP, antibody 
positive swine sera showed characteristic cellular surface 
staining patterns (Figure 1C), which were indistinguish- 
able from those of REBOV-infected cynomolgus monkey 
sera and REBOV-GP immunized rabbit sera (data not 
shown). Antibody-negative swine sera showed no reac- 
tion (Figure ID). In Bulacan, 104 (71.2%) and 115 
(78.8%) of the 146 swine sera showed positive reactions 
in the NP- and GP-specific IFA, respectively. In Pangasi- 
nan, each 54 (78.3%) of the 69 swine sera showed posi- 
tive reactions in the NP- and GP-specific IFA. In total, 
158 (73.5%) and 169 (78.6%) of the 215 swine sera col- 
lected at the affected farms were REBOV-NP and -GP 
antibody positive in the IFA, respectively (Table 1). On 
the other hand, none of the 98 swine sera collected in 
Tarlac in the Philippines showed any positive reaction 
(Table 1). 



Detection of REBOV neutralizing antibodies in swine 

The VSV-pseudotype bearing REBOV-GP efficiently 
infected Vero E6 cells which are known to be susceptible 
to REBOV infection. The infectious titer of the VSV- 
pseudotype reached 3.6 x 10 6 IU/mL when measured on 
Vero E6 cells. Infection with the VSV-pseudotype was 
neutralized with rabbit serum to REBOV-GP (data not 
shown). Swine sera were analyzed in the NT using the 
VSV-pseudotype. In the NT, none of the 49 swine sera 
collected in Japan neutralized the infection with the 
VSV-pseudotype on Vero E6 cells at serum dilutions of 
1 in 100. These were therefore considered to be REBOV 
NT antibody negative. In the affected farms, 108 (74.0%) 
of the 146 swine sera in Bulacan and 46 (66.7%) of the 
69 swine sera in Pangasinan showed positive reactions in 
the NT, respectively at a serum dilution of 1 in 100. In 
total, 154 (71.6%) of the 215 swine sera in the affected 
farms showed NT antibody positive (Table 1). NT titers 
of the sera were then obtained for the 34 NT antibody 
positive sera in the affected farms. NT titers ranged be- 
tween 100 and 12,800 with average of 790 and median 
of 400 (data not shown). In Tarlac, none of the 98 swine 
sera showed any positive reaction (Table 1). 

Detection of REBOV-NP and -GP specific antibodies in 
swine using IgG-ELISA 

IgG antibodies to the REBOV-NP and -GP were 
detected by recombinant REBOV protein-based IgG- 
ELISA. The OD values to the NP and GP of each 
serum sample at a dilution of 1 in 100 were plotted 
(Figure 2A and B), and ROC and TG-ROC curves were 
also drawn (Figure 2C to F). Since the OD values of 
negative samples were very low, the cut off values 
defined as the values at intersection points of the sensi- 
tivity curve and specificity curve were 0.077 and 0.104 
for REBOV-NP and -GP, respectively. At the cut off 
values, sensitivity and specificity of the assay are over 
95% in both IgG-ELISAs (Figure 2E and F). In Bulacan, 
115 (78.8%) and 119 (81.5%) of the 146 swine sera 
showed positive reactions in the NP- and GP-specific 
IgG-ELISA, respectively. In Pangasinan, 62 (89.9%) and 
46 (66.7%) of the 69 swine sera showed positive reac- 
tions in the NP- and GP-specific IgG-ELISA, respect- 
ively. In total, 177 (82.3%) and 165 (76.7%) of the 215 
swine sera collected at the affected farms were 
REBOV-NP and -GP antibody positive in the IgG- 
ELISA, respectively (Table 1). Conversely, 1 of the 98 
swine sera collected in Tarlac in the Philippines showed 
a positive reaction in the GP specific IgG-ELISA and 1 
of the 49 swine sera collected in Japan showed a posi- 
tive reaction in the NP- and GP-specific IgG-ELISA 
(Table 1). However, these results are considered to be 
false positives because the OD values were close to the 
cut off level defined by the assay, additionally all the 
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Figure 1 Detection of REBOV antibodies in IFA. Immunofluorescence staining pattern of HeLa cells expressing REBOV-NP (A) and REBOV-GP 
(C) with a REBOV antibody positive swine serum at a dilution of 1 in 160. Negative staining pattern of HeLa cells expressing REBOV-NP (B) and 
REBOV-GP (D) with a REBOV antibody negative swine serum (B and D). 
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Table 1 Detection of REBOV-antibodies in swine 





NP 




GP 




NT 
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IgG-ELISA 


IFA 


Bulacan 


115/146 (78.8%) 


104/146 (71.2%) 


119/146 (81.5%) 


115/146 (78.8%) 


108/146 (74.0%) 


Pangasinan 


62/69 (89.9%) 


54/69 (78.3%) 


46/69 (66.7%) 


54/69 (78.3%) 


46/69 (66.7%) 


Tarlac 


0/98 (0%) 


0/98 (0%) 


1/98 (1.0%) 


0/98 (0%) 


0/98 (0%) 


Japan 


1/49 (2.0%) 


0/49 (0%) 


1/49 (2.0%) 


0/49 (0%) 


0/49 (0%) 



samples collected in Tarlac and Japan tested using the 
alternative serological assays were negative. 

Agreement between the serological assays detecting 
REBOV antibodies 

We next examined the agreement between the sero- 
logical assays used in the present study. The data 
obtained in each assay were binarized to either positive 
or negative result. These binarized data were then statis- 
tically analyzed by pair-wise comparisons in nonpara- 
metric one-way ANOVA. This analysis showed that 
there were no significant differences between the sero- 
logical assays (p > 0.05, data not shown). The result sup- 
ports the contention that all of the five serological assays 
used in the present study posses similar discriminatory 
capacity. Thus we are confident that many REBOV anti- 
body positive samples showed a positive reaction in all 
of the five assays. 



Discussion 

ZEBOV and Marburg virus were identified in African 
fruit bat species [19,20]. More recently, Marburg virus 
was successfully isolated from a fruit bat, Rousettus 
aegyptiacus [21], suggesting that fruit bats are reservoir 
animals of filoviruses. In the Philippines, we have re- 
cently demonstrated that an Asian fruit bat, Rousettus 
amplexicaudatus, has antibodies to REBOV [22]. Since 
the REBOV genome has not yet been detected in the 
bat, conclusive evidence that the bat species is a reser- 
voir, or one of the reservoir animals, of REBOV is not 
available. Nevertheless, it is possible that REBOV was 
transmitted to swine from these bats since these bats in- 
habit many areas of the country, including the regions 
around the affected facilities both in Pangasinan and 
Bulacan. 

In this study, we have aimed to clarify how REBOV in- 
fection was spread among swine during the REBOV 
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Figure 2 Distribution of OD values and ROC and TG-ROC analyses of the IgG-ELISA specific for REBOV-NP and -GP. OD values to the NP 
(A) and GP (B) of each serum sample at a dilution of 1 in 100 were plotted for swine group showing either IF and NT antibody positive or 
negative. ROC and two graph-ROC (TG-ROC) curves were analyzed and drawn using Stat Flex software. 
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epizootics. We have tested 215 swine sera collected from 
the REBOV affected farms using IFA and IgG-ELISA 
specific for REBOV-NP and -GP, and NT. NT is the gold 
standard of serological assay in many virus infections. 
However, since REBOV needs to be cultured in high- 
containment laboratories, we performed an alternative 
NT using the VSV-pseudotype bearing REBOV-GP. This 
avoided the use of infectious REBOV, enabling the work 
to be carried out at low containment. Previously it has 
been shown that VSV-pseudotype bearing ebolavirus GP 
mimicks ebolavirus infection [23], Approxymately 70% 
of the swine sera from REBOV affected farms were 
REBOV antibody positive. This indicated that swine are 
susceptible to REBOV infection. Unfortunately, we could 
not analyze the IgM antibody responses in swine, since 
after the sera were heat inactivated the gamma globulin 
fractions were precipitated with ammonium sulfate, and 
reconstituted in PBS prior to be testing. An indication of 
the IgM responses to REBOV would have provided evi- 
dence of a recent infection. Thus, it is still unclear if 
REBOV infection was spread during epizootics or 
whether a population of the animals in the farms was 
infected with REBOV prior to the epizootics. The swine 
not associated with the epizootics, in Tarlac, are consid- 
ered to be free from REBOV infection. These samples 
were collected in 2010, over 2 years after the epizootic, 
and from animals born after the epizootic. Moreover, we 
could not analyze the swine specimens near the affected 
farms in 2008. Thus, in this study, it is not clear if 
REBOV infection in 2008 was limited in the affected 
farms. Further study is necessary to conclude if the 
swine population in the Philippines is generally free from 
REBOV infection. 

Recently, it has been shown that the experimental infec- 
tion of swine with REBOV alone resulted in subclinical in- 
fection with rapid clearance of the virus [24]. Alternatively, 
ZEBOV has been shown to replicate to high titers in ex- 
perimentally infected swine and to cause severe lung path- 
ology resulting in transmission of the virus to naive animals 
[25]. Thus, swine has been shown experimentally to be 
highly susceptible to ZEBOV infection. Furthermore, some 
amino acid mutations in NP and/or VP24 in ZEBOV 
resulted in adaptation of the virus to guinea pigs and mice 
[26,27]. Thus, we cannot rule out the possibility that muta- 
tions introduced in the REBOV genome during serial trans- 
mission in swine will result in adaptation of the virus to 
swine in future. In this regard, a regular serological survey 
of REBOV infection in swine in the Philippines is desirable. 
The serodiagnosis systems presented in this study might be 
useful for such a survey. 

Conclusions 

The high prevalence of REBOV infection in swine at the 
affected farms in 2008 suggests that swine are susceptible 



for REBOV infection. The multiple serological assays used 
in the study are thought to be useful for future surveil- 
lance of REOBV infection in swine in the Philippines. 

Methods 

Swine serum specimens 

A total of 215 swine sera were collected from two 
REBOV affected pig farms, located in Pangasinan and 
Bulacan in 2008 (Figure 3). Of these, 146 sera were col- 
lected from swine at the farm in Bulacan, and 69 were 
collected from those in Pangasinan. Swine samples in 
the affected farms were collected under quarantine of 
the Philippines. The sera were kept frozen at the 
Research Institute for Tropical Medicine (RITM) in the 
Philippines until use. Ninety-eight swine sera were col- 
lected from July to September 2010 from swine aged 
between 2 and 20 months (median of 4.5 months) in 
Tarlac in the Philippines, where no swine epizootic has 
been documented. The swine specimens at Tarlac were 
collected and used under approval of IRB (No. 2009-018) 
of RITM, and informed consent was obtained from the 
farm owners. Forty-nine swine sera collected at slaugh- 
terhouse in 2006-7 in Japan, kindly supplied by Dr. T-C 
Li at the National Institute of Infectious Diseases, were 
used as REBOV antibody negative control sera. Swine 
sera were inactivated at 56°C for 30 minutes, and then 
the gamma globulin fractions were precipitated with 
ammonium sulfate and reconstituted in phosphate buf- 
fered saline (PBS). 

Immunofluorescence assay (IFA) specific to REBOV-NP 
and -GP using HeLa cells expressing the recombinant 
protein 

The entire cDNA of REBOV-GP ORF was amplified 
from the swine lymph node specimen used for amplifica- 
tion of the cDNA of NP by RT-PCR using the primers 
REBOV-GP/F (5'-CGA AGC TTC GAA CAT GGG 
GTC AGG ATA TCA ACT-3') and REBOV-GP/R (5'- 
CGA AGC TTC AAC ACA AAA TCT TAC ATA TAC 
AAA G-3') (the Hind III site is underlined). The 
complete nucleotide sequence of the amplicon was 
determined to be identical to that of Reston 08A (Gen- 
Bank accession number FJ621583), indicating that the 
specimen was collected from the same farm as that of 
sample group A from which Reston 08A was isolated 
[17]. The Hind III fragment of the amplicon was subse- 
quently cloned into pKS336 [28], and the pKS336 plas- 
mid with REBOV-GP, pKS336-pREBOV-GP, was used 
for transfection. HeLa cells were transfected with 
pKS336-pREBOV-GP using a FuGENE HD (Roche Diag- 
nostics, Mannheim, Germany) according to the manu- 
facturers instructions. The cells transfected with the 
plasmids were selected with 3 ug/mL of blasticidin S- 
hydrochloride (Invitrogen, CA, USA) in Dulbeccos 
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Figure 3 A topographical map showing locations in the Philippines where swine serum samples were obtained. Swine serum samples 

were collected at farms in three regions in the Philippines, Bulacan, Pangasinan and Tarlac. These regions were showed the numbers of 

specimens and percentage of positivity using NT. The three locations and Laguna, where the Ferlite monkey facility used to be located, are 

shown on a Google map (©201 1 Google Map data ©201 1 AND Europa Technologies), 
v J 



modified Eagle medium (DMEM) supplemented with 5% 
fetal calf serum (FCS). The selected cells were analyzed 
for the expression of REBOV-GP in IFA using anti- 
REBOV-GP rabbit serum. The cells expressing 
REBOV-GP were spotted on multiwell slide glasses, 
fixed in acetone, and used as antigens for IFA specific to 
REBOV-GP. The IFA specific to REBOV-NP was per- 
formed as described previously [18], with slight modifi- 
cations. Mock-transfected HeLa cells were used as 



negative control antigens in the IFA. The antigen slides 
were incubated with serially diluted swine sera under 
humidified conditions at 37°C for 1 hour. The antigen 
slides were washed in PBS, then reacted with rabbit anti- 
pig IgG conjugated with fluorescein isothiocianate 
(FITC) (Bethyl, TX, USA) at a dilution of 1 in 100 at 37° 
C for 1 hour. The slides were washed in PBS, covered 
with cover glasses, and examined for staining patterns 
under a fluorescence microscope fitted with appropriate 
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barrier and excitation filters for FITC visualization. The 
antibody titer in the IFA was determined as the recipro- 
cal of the highest dilution showing positive staining. 

Generation of vesicular stomatitis virus (VSV)-pseudotype 
bearing REBOV-GP 

Generation of VSV-pseudotype bearing REBOV-GP was 
performed as described previously [29-31]. Briefly, 293 T 
cells were transfected with pKS336-pREBOV-GP and 
cultured for 24 hours, and then the cells were infected 
with VSVAG* (kindly supplied by Prof. M.A. Whitt, Uni- 
versity of Tennessee) [30,31] and cultured for 24 hours. 
The culture supernatants were then collected, filtered 
through a 0.22 (im-pore-size filter, and stored at -80°C 
until use. The infectious titer (infectious unit, IU) of 
VSV-pseudotype on Vero E6 cells was determined by 
counting the number of Green Fluorescent Protein 
(GFP)-expressing cells. 

Neutralization test (NT) 

The sera were diluted twofold from 1 in 100 with 
DMEM containing 5% FCS and 1,000 IU of VSV- 
pseudotype bearing REBOV-GP. The mixture was incu- 
bated for 1 hour at 37°C, then inoculated onto Vero E6 
cells seeded on 96-well plates. The numbers of VSV- 
pseudotype infected cells were determined by counting 
of the number of GFP-positive cells according to the 
methods described previously [23,29,30]. NT titers of 
the tested sera were defined as the reciprocals of the 
highest dilutions at which more than 50% inhibition of 
infectivity was observed. Serum samples were considered 
to be NT antibody negative when less than 50% inhib- 
ition of infectivity was observed at a dilution of 1 in 100. 

Recombinant baculoviruses expressing recombinant 
REBOV proteins 

cDNA of REBOV nucleoprotein (NP) open reading frame 
(ORF) from a swine lymph node specimen collected at the 
farm in Bulacan was amplified by reverse transcription 
polymerase chain reaction (RT-PCR), and the nucleotide 
sequence of the amplicon was determined to be identical to 
that of Reston 08A. The cDNA was subcloned into pGEM- 
Teasy (Promega, WI, USA) and used for the following PCR 
template. PCR was performed using the plasmid clone to 
add BamHl linker with primers REBOV-NP/F (5'-GGG 
CTA GCG GAT CCA AGT CGA TAT GGA TCG TGG 
GAC C-3') and REBOV-NP/R (5'-TTG CGG CCG CGG 
ATC CCT GAT GGT GCT GCA AGA TTG-3') (the re- 
striction site is underlined). A itamHI-digested fragment of 
the plasmid was then subcloned into pAcYMl-C-His plas- 
mid, a derivative of pAcYMl plasmid [32] carrying eight 
histidine coding sequences just downstream of the BamHl 
site, to construct pAcYMl-His-pREBOV-NP. A recombin- 
ant baculovirus, Ac-His-pREBOV-NP, was generated using 



a previously described method [33]. The recombinant 
baculovirus, which expresses the ectodomain of REBOV- 
glycoprotein (GP) with a histidine-tag at its carboxyl 
terminus [22], was used to prepare recombinant REBOV- 
GP for IgG-ELISA. A baculovirus (Ac-AP) that lacks the 
polyhedorin gene was used to prepare a negative-control 
antigen in insect cells. 

Expression and purification of recombinant NP and GP of 
REBOV in baculovirus 

Tn5 insect cells were infected with Ac-His-pREBOV-NP 
and Ac-His-REBOV-GP and then incubated at 26°C for 
72 hours and 48 hours, respectively. The cells were 
washed in PBS, lyzed in PBS containing 1% NP40 and 
8 M urea, and then clarified by centrifugation at 
8,000 rpm for 10 min. The supernatant fraction was col- 
lected, and recombinant REBOV-NP and REBOV-GP 
were purified using a Ni2 + -resin purification system 
(QIAGEN, Diisseldorf, Germany) according to the man- 
ufacturers instructions. The expression and purification 
of REBOV-NP and REBOV-GP were confirmed in 10% 
sodium dodecyl sulfate polyacrylamide gel electrophor- 
esis (SDS-PAGE) and Coomassie Brilliant Blue R-250 
staining Tn5 cells infected with Ac-AP were processed 
similarly and used as negative control antigen in the 
ELISAs. The purified antigens and negative control anti- 
gen were kept at -80°C until use. 

IgG-ELISA 

Swine sera were analyzed for the presence of antibodies 
to REBOV-NP and REBOV-GP in IgG-ELISA, essentially 
as described previously [22,34,35]. Briefly, half of the 
wells of the ELISA plates (Falcon, NJ, USA) were coated 
with a predetermined optimal quantity (approximately 
100 ng/well) of each purified recombinant protein, and 
the rest of the wells were coated with the negative con- 
trol antigens. After washing the plates three times in 
PBS containing 0.05% Tween-20 (T-PBS; SIGMA- 
ALDRICH, MO, USA), each well was incubated with 
200 \iL of T-PBS containing 5% skim milk (MT-PBS; 
Yukijirushi, Hokkaido, Japan) for 1 hour at 37°C. After 
washing the plates three times in T-PBS, the antigen- 
coated and negative-control-antigen-coated wells were 
inoculated with the test samples (100 (iL/well) in MT- 
PBS at a dilution of 1 in 100, incubated for 1 hour at 37° 
C, and washed in T-PBS. Then, each well was incubated 
with a Protein A/G conjugated with horseradish perox- 
idase at a dilution of 1 in 1,500 (PIERCE, IL, USA) in 
MT-PBS for 1 hour at 37°C. The plates were washed 
three times, and 100 \xL of ABTS solution (Roche Diag- 
nostics, Mannheim, Germany) was added to each well. 
The plates were incubated for 30 min, and the optical 
density (OD) was measured at 405 nm with a reference 
at 490 nm. The adjusted ODs for each tested sample 



Sayama et a I. BMC Veterinary Research 2012, 8:82 
http://www.biomedcentral.eom/1 746-61 48/8/82 



Page 8 of 9 



were calculated by subtracting the ODs of the negative- 
control-antigen-coated wells from those of the corre- 
sponding REBOV- antigen- coated wells. 

Statistical methods 

Sensitivity, specificity and predictive values for positive 
and negative tests were calculated by standard methods. 
Sensitivity and specificity were defined as the probability 
that the target assay result was positive when the IFAs 
specific to REBOV-NP and -GP and NT showed positive 
and the probability that the target assay result was nega- 
tive when the IFAs and NT showed negative, respect- 
ively. Receiver operating characteristics (ROC) and two 
graph-ROC (TG-ROC) curves were analyzed using Stat 
Flex software (Artech Co. Ltd., Osaka, Japan) [36,37]. To 
examine the agreement between the serological assays 
used in the present study, we employed Friedman test 
with Dunn's post-hoc test [38,39] using GraphPad Prism 
using pair-wise comparisons in nonparametric one-way 
ANOVA. 
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